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P3 . Abstract 

^ I We analyze here the energy states and associated wave functions available to a particle acted 

upon by a delta function potential of arbitrary strength and sign and fixed anywhere within a 

^ ' one-dimensional infinite well. We consider how the allowed energies vary with the well's width and 

^. 

ff^ , with the location of the delta function within it. The model subtly distinguishes between whether 



Cn ' the delta function is located at rational or irrational fractions of the well's width: in the former 

o 
o 



C^ ' case all possible energy eigenvalues are solutions to a straightforward dispersion relation, but in 



the later case, to make up a complete set these 'ordinary' solutions must be augmented by the 
addition of 'nodal' states which vanish at the delta function and so do not 'see' it. Thus, although 

X" 

^ ■ the model is a simple one, due to its singular nature it needs a little careful analysis. The model, 

of course, can be thought of as a limit of more physical smooth potentials which, though readily 
succumbing to straightforward numerical computation, would give little generic information. 
PACS numbers: 03.65.-w, 73.21.Fg, 01.40.-d 



I. INTRODUCTION 

Calculating the energy states for the motion of a particle in a 'simple' one dimensional 
infinite square well finds its way into all standard quantum mechanics texts. But because 
the potential representing the walls is sharp and infinite the model has subtleties. For 
instance a consistent definition of a momentum operator for it is challengingii^. Another 
intriguing property of the infinite well is that for certain initial states the time development 
of the spatial probability density can show fractal behavior=2. There is also the topic of wave 
packet revival in such a well^'^. And the time dependence of a particle's wave function in 
a suddenly expanded well can show interesting features^. Clearly, when sharp edges and 
infinities are involved, care is important. 

In this paper we compound the singularity by adding within the well a fixed delta function 
potential of arbitrary sign, strength and location. For this model, due to the containing 
walls, all solutions to the Schrodinger equation are bound states, but the model does have 
its niceties. We find for instance that although, as expected for the bound states of any 
one-dimensional system, all of its energy eigen-states are non-degenerate, should the delta 
function be located at a rational fraction of the well's width, then a subset of these energies 
approach in the limit of strong attraction or repulsion, double degeneracy. 

In addition to examining several other properties and limits of the model we have 
computed — for both attraction and repulsion — the ground state energy of the system as a 
function of the location of the delta function. For the case of attraction this energy has a 
minimum when the force center is located at the center of the well. For a wide box this 
minimum is nearly fiat, but it becomes sharper when the well width is reduced to sizes of 
order of the spatial decay length characterizing the 'molecular' size of the bound state of an 
attractive delta function in free space. 

The model has not, of course, escaped attention. PatiF and Atkinson and Crater* give 
analyses for the case that the delta function potential is located precisely at the center of the 
well. Bera and his co-workers^ place it anywhere within the well and utilize perturbation 
theory. Joglekar— computes the energy eigenvalues when it is located at irrational multiples 
of the well's width and considers weak and strong coupling limits. Epsteinii uses the well 
with the delta function located at the center in order to gain insight into the interesting 
question posed by Wigner— as to whether the energy levels of a hydrogen atom should be 



derivable from second order perturbation theory in view of the fact that it is known to be 
proportional to the fourth power of the electron charge and therefore to the square of the 
Coulomb interaction. 

We were motivated to consider the model by the long standing problem of the effects of 
containment on molecule states and energies dating back to work by Michels^^,^ Sommerfeld^ 
and their co-workers, in the early heyday of quantum mechanics, applying approximations 
to study what might be said to be the canonical problem — that of a hydrogen atom with 
its nucleus fixed at the center of an impenetrable sphere. Other more recent papers discuss 
this model with various methods of approximation. See for instance references'^"— and, 
especially22 which includes a short review of work on this model up to 1987. More recently 
Changa and co-workers2i have applied perturbation theory to estimate the energies when 
the hydrogen nucleus is shifted from the center of the the spherical box. 

The particular advantage of the model we consider here is that no approximations are 
necessary to discuss the energy levels and their dependence either on the well width or on 
the location of the force center. It would be gratifying to be able to do the same with other 
model interactions. For instance one might try a 'one dimensional Coulomb' repulsion or 
attraction, v{x) ~ l/|a^|, a; 7^ 0. But it is a sobering fact that even with no walls present, so 
that X is any real number bar zero, there seems to be no final agreement for that model as 
to the available energy levels and corresponding eigenstates^Si"^— . 

The paper is arranged as follows. We give the general solution in Sec. [TTl Its nature 
depends upon whether the ratio i/L is a rational or irrational number, where i is the delta 
function's position within a well of width L. When that ratio is irrational all solutions 
are straightforward 'ordinary' solutions, but when it is rational a complete set of solutions 
includes both the ordinary ones together with those we choose to call 'nodal' solutions. For 
both cases we confirm in Sec. IIIII mutual orthogonality between all solutions corresponding 
to different energies. In Sees. [IV] and |V] we consider the limits of weak and strong coupling 
for any location of the delta function potential within the well. In taking these limits, when 
i/L is rational care must be exercised to include both the nodal and ordinary solutions. In 
Sec. I VII we discuss solutions for the special symmetric case that the delta function is located 
precisely in the middle of the well, for which i/L is the rational number 1/2. In Sec. I VI 1 1 we 
consider in some detail the ground state of the model, especially with respect to its energy 
dependence on the relative location i/L, for various signs and strengths of the delta function 



potential. This is plotted in Fig. |2l 

II. SOLUTION 

A particle of mass m free to move within a one-dimensional well with sides at a; = and 
X = L has energy eigenstates 

<W = y|sm(^), iV=l,2,3,..., (1) 

and corresponding energies 

For odd values of N these states are symmetric about the middle (x = L/2) while those for 
even A^ are antisymmetric. 

A particle constrained by no walls but acted upon by a one-dimensional delta function 
potential located ai x = i, namely 

v(x) = -X6{x-e), (3) 

has the associated Schrodinger equation 

-— $"(x)-A(5(x-£)^(x) = E^(x). (4) 

where, if \E'_ and \E'+ are solutions to the left {x < i) and right {x > £) of the delta function, 
the conditions to be satisfied at the potential are 

(^+(^) - ^-(^)) = A^_m and ^_m = ^+(i). (5) 

For an attractive interaction (A > 0) the single bound state is represented by 



where 

A. 4 m 

mX 
is a measure of the 'size' of the bound state and the energy is 

mX' _ h' 1 



Whichever the sign of A, solutions to @ also include a continuum of positive energy unbound 
states, the scattering states. 

Now consider putting the potential (^ inside a box of finite side L. Then all states of 
the system are bound under the combined action of the walls and delta function. For this 
model, solutions must satisfy Eq. (jl]), must vanish at the walls and must obey conditions 
([5]). For X ^ i solutions of Eqs. (|3]) and (jl]) are linear combinations of exp{±ikx), with 
energy E = Ti k"^ /2m. To be continuous at x = £ and vanish at x = and x = L solutions 
will have the form 

{ip+(x) = C sinikiL — x)) sm(ki) ior i < x < L 
(9) 
^_ (x) = C sm{kx) sm{k{L - i)) for < x < £, 

where C is a constant which may depend on k, L and i. 

Applying the first of boundary conditions ([5]) gives the energy dispersion relation. It is 

kA sm{kL) = 2 sin(H) sin(A;(L - £)), (10) 

where A is given by ([7]). One set of solutions to flTOj) is given by 

sin(A;L) = 0, sin(H) = 0, sin(fc(L - £)) = 0, (11) 

where any one of these equations implies the other two. Here k must meet the two-fold 

condition 

k = -- = ^ so that - = -, 12 

L £ L n 

where n and p are positive integers. We shall assume that the ratio p/n has been reduced 
to its primitive form — that is to say with no common factors other than unity — so that, 
because i/L = p/n = pj/{nj), we can write all of these nodal solutions ^Ij^"-\ and their 
energies, for any such given value of i/L as 

-sin( 
with wave numbers given by 



V'S(^) = \/;^sin(-^j and E.,, = — (^) , j = l,2,3..., (13) 



kL = jn-K. (14) 

The particular feature of these solutions is that they have nodes at the location of the delta 
function potential. They are that subset of the standing wave solutions ^ for a free particle 



in a well with this property. Because they vanish at x = i they satisfy both conditions iQ. 
As an example, when i/L equals 1/3 (or 2/3) the states are ^3"- (a;) = \/2/L sm{3njx/L) 
with energies E^j = {Ti /2m)(37rj'/L)^ where j = 1, 2, 3 . . .. 

When s\.n{kL) does not vanish Eq. ( TTOl) can be re-expressed as 



^^^^^^^^ .^^-e, ^ (-) 



or, equivalently, as 



— — = tan — - 2 }-^-^ where u = 2 1 (16) 

2 \L J V 2 / sin(A;L) ^ L ^ ^ 

and i varies from to L as /i ranges from —1 to 1. We choose to call solutions for k to either 
of these versions of the dispersion relation (together with their associated wave functions) 
ordinary solutions. 

Should i/L be a rational number, complete information about the nodal states and their 
energies (all positive) are given by Eqs. (1131) and (TT^ . The ordinary states exist for all 
values of £ between and L, with energies E = h'^k'^ /{2m). For positive energies the allowed 
values of k are real-valued. In that case, solutions for kL to Eq. (IT5l) occur at intersections 
with the expression on the right-hand side of a straight line passing through the origin with 
slope (A/L). Similar comments apply to Eq. flT6|) in terms of kL/2 and slope 2A/L. For 
negative energies, should there be any, the wave number is fc = Ik, where k is real. In that 
case the energy is E = —h k? /{2m) and, letting A; — )• i/t in the pair ( TT5l) and ( TT6l) and in (I9l), 
one must solve either of the equations 

with corresponding wave functions 

i){x) = { (18) 

D sinh(Ka;) smh.{K{L ~ i)) foi < x < i . 

III. ORTHOGONALITY 

We choose to confirm by direct calculation that for this one-dimensional problem all 
solutions corresponding to different energies are, as expected, mutually orthogonal. 



Note first that, because the nodal solutions are free particle states in a simple well, it 
follows at once that those for different energies are mutually orthogonal. Then too these 
nodal solutions, Eq. ( TT^ . are orthogonal to the ordinary states. To see this, write the nodal 
states in condensed notation as 



iIj^'^'{x) = y — sin(z/x), where z/ = — — = — — . (19) 

\ J-/ -L/ -t 

Now, for k ^ and k ^ p, 

usin{k£) 



sin(i/x) sm{kx) dx 



\pj 



k"^ — v"^ 
and 



/ sm(z/x) sm(A;(L — xj) dx = — (— j'^"' — . 

Jn k^ — v^ 

Then using definition (|9]) gives 

V'(x)V^(")(x)dx= \ ip_{x)ip^''\x) dx + f ij+{x)ij^''\x) dx = 0. (20) 

Jo Je 

Orthogonality also holds between any t/^^"^ and any ordinary solution ip with negative energies 

E = —Ti^H? /{2m), where k is a solution to equation ( I17p . 

Finally, the ordinary solutions (|9]), where k satisfies the dispersion relation (TT5|) . are 

orthogonal to each other for different values of k. For instance, including subscripts to 

indicate the /c-dependence of wave functions, ipk,±{x) and normalization constants C^, we 

find by direct calculation, that when k ^ k', 

L f-i i-L 

i)k{x)i)k'{x)dx = / ilJk-{x)ipk'-{x)dx+ / ilJk,+ {^)'ipk',+ {^)dx 

Jo Je. 

_ CkCk' 
~ 2 

where, doing the integrals and rearranging 



(A-B) 



A = —^—-\sm(k(L-i))sm(k'(L-i))sm((k-k')i) 
[k — k') '- 

+ sin(fc£) sin(A;'£) sin((A; - k'){L - £))] 

and B = —^^--\sm(k(L-i))sm(k'(L-i))sm((k + k')i) 
[k + k') '- 



+ sin(H) sin(A;7) sin((A; + k'){L - 
Then, by using the eigenvalue equation (IT5|) and trigonometric identities one finds that 

A = - sm{kL) sm{k'L) = B, 



so that, when k ^ k', 

Mx)^k'ix)dx = 0. (21) 



IV. WEAK COUPLING 

The need to include both the nodal and ordinary solutions to make up a complete set 
can be seen by considering the limit of weak coupling for which |A| is small and |A|/L large. 
In that limit one would expect, whatever the value of ^ G (0, L), that solutions ([9]) and their 
energies would approach those of a free particle in a well, Eqs. ([1]) and ([2]), no more and no 
less. 

Consider solutions with positive energies. For arbitrary i G (0,L), in addition to the 
nodal states — which exist only when i/L is rational — we have the ordinary states. For them 
kL/n is not an integer and the dispersion relation can be written, say, as f lTSj) . Solutions 
for kL to either equation are intersections of a straight line passing through the origin with 
slope A/L, large positive for attraction or large negative for repulsion. These intersections 
occur only near any divergences of the right-hand side (RHS) of Eq. f llSp . i.e. near the 
zeroes of sin(A;L). Then, setting kL = Ntt + e, where e is small but of either sign, shows 
that near these zeroes we have, approximately, 

-(2/e) sm\NTTi/L) + ^ a„(A^) e" (for A^ = 1, 2, 3 . . .) 
RHS ~ { "^° (22) 

(e/2)(l - /i2) + b{e) (for N = 0, e < tt), 

where 6(e) is monotonically increasing as a function of e in the interval [0, vr) and /i = 
2{i/L) - 1. 

Suppose first that i/L is irrational. In that case the function sin(A^7r£/L) doesn't vanish 
for any non-zero integer A^. Then, using the first of these two equations in flTSj) . with 
kL = N7T + e, gives (for A^ = 1, 2, 3, . . .), 

sm\N7re/L) , , ...,.x 

e = —2 -— - — h terms m higher mverse powers of (A/L). 

N'k{A/L) \ I j 

Using the second of equations ( 122|) in ( TT5|) verifies that there is no solution in this limit 
for A^ = 0. To summarize: when i/L is irrational, solutions for kL to ( ITSll — or, what is 
equivalent, ( IT6|) — will, for ever weaker coupling (|A|/L — ;■ oo), approach A^vr for attraction 
or repulsion, where A^ = 1,2,3..., so making up the full set of energies for a free particle 
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in a well. Also, it is easy to see, in the limit kL = Ntt + e with e — )■ 0, that solutions ([9]) 
approach those of the free particle in a well, Eq. ([1]). 

If instead i/L should take the rational value p/n, which we assume has been reduced by 
eliminating all common factors, then sin(A^7r£/L) will vanish when and only when N = j n, 
j = 1, 2, 3 . . .. There are in the limit, therefore, no ordinary solutions for these values of N, 
but we do have the nodal states, Eqs. ( !T3|) and (fT^ . In addition to these nodal states, as we 
know from the argument just given for irrational i/L, for large |A|/L there exist ordinary 
solutions for all other integral values of A^ again making up the full set of solutions for a free 
particle in a one-dimensional well. 

For negative energies, similar arguments show that when |A|/L is large there are no 
solutions to either of Eqs. flT7|l . 

V. STRONG COUPLING 

Now consider, for any i G (0,L), the strong coupling limit, for which the magnitude of 
A/L is small. Graphical argument shows that in this case, and only for attraction (A > 0), 
there is a single state of negative energy. This is easily verified by considering Eq. f fT7|) 
for small positive A/L. For any value of /i between —1 and +1 the right-hand side of that 
equation is positive and monotonically increasing as a function of positive kL, from zero 
asymptotically to the value +1. Thus, the single solution approaches the value kA = 1, 
giving the limiting energy value E = —h h? /{2m) = —Eb- 

Except for the single bound state when the potential is attractive, for strong coupling the 
behavior of the oscillatory positive energy solutions to Eq. ^ is controlled to some extent 
by the product Xip{i) in the potential energy term, for which one expects that whatever its 
sign, as the magnitude of A increases any ordinary solution evaluated at the potential, ipi^i), 
will decrease in compensation. Of course when i/L is rational, the corresponding nodal 
solutions strictly vanish at the delta function and so are independent of the magnitude of 
A, so long as it is non-zero. 

For the positive energies, when \A\/L is small, we can say, from equation ( 1T0|) . that 
sm{ki) sm{k{L — i)) is also small, so that, whatever the sign of A, 

either (a) ki = riiir + Si or (b) k{L — i) = n2n + 62, (23) 



where rii and n2 are any positive integers and ei and £2 are small and vanish as |A|/L — )■ 0. 
Note that, in the limits of vanishing ei and £2, cases (a) and (b) cannot both hold, for then, 
according to OH]), the wave function would strictly vanish. 

Suppose first that d/L is irrational. Then the two cases given in ( 12^ generate, in the 
limit, distinct energies and do not overlap for any values of ni and 722: If this were not so, 
then for some integers rii and 722 we would have kL = niir and k{L — i) = n2TT for the same 
value of k so that i/L would be rational. 

Cases (a) and (b) of (!23|) give, in the limit of vanishing A/L, nothing other than the 
energy levels of a free particle in separate wells of width i and {L — €). We can underline 
this by deriving the wave functions in the strong coupling limit. For case (a), using in Eq. 
on]) the first of expressions fl25]) and taking the limit Si — )■ 0, gives 

.0 for i< x< L 

^ sm{ni7rx/i) for < x < i. 



And, similarly, for case (b): in the limit £2—^0 these states approach the ordinary well states 
constrained to the interval [i, L] with quantum number ^2. That the energies approach those 
of free particle motion in separate wells of widths i and L — £ as |A|/L — ;■ (when i/L is 
an irrational number) has been pointed out by Joglekar— . 

There are subtleties when i/L is rational, with i/L = p/n, where we assume as usual 
that all common factors have been eliminated from numerator and denominator. In that 
case, for all values of A, we have the nodal solutions, Eqs. (fT3|) and ([HD, with wave numbers 
given by kL = jmr, where j = 1, 2, 3 ■ ■ ■ . And, as before, in addition to those we have the 
ordinary solutions for which sin(fc^) sin(A;(L — i)) is small. For these, from Eq. (!23ll . 

either (a) kL = riin (n/p) + Ea or (b) kL = n2n n/{n — p) + Eb, (24) 

where Ea and Eb are small. In this case possibilities (a) and (b) are both satisfied whenever 
n2/ni = {n~p)/p . But that instance gives — as we have said — a vanishing wave function and 
so must be ruled out as a solution. Another property is that when rii = jp or ^2 = {n —p)j, 
where j = 1, 2, 3 ■ ■ ■ , they are converging to nodal state energies. But note that, however 
small |A|/L may be, so long as it is non-zero these ordinary states are distinct from nodal 
states and orthogonal to them. In summary: when i/L = p/n (reduced) then the set of 
limiting values for kL as |A|/L approaches zero comprises the nodal values together with all 
distinct (non-duplicating) limits in Eq. (IMj) . 
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As an example we have considered the case C./L = 2/5. Its nodal states have wave vectors 
given by kL/ir = 5j, j = 1, 2, 3, . . .. In addition to these we have the ordinary states. In the 
limit |A|/L — ;■ these are given by the zeroes of the right-hand side (RHS) of Eq. (TT^ (and 
fITB]) ) with i/L = 2/5. This is shown as the dashed lines in Fig. [T] as a function of kL/ir 
for values up to kL/n = 9. The solid curves show the same function plotted for the nearby 
value i/L = 0.415 = 83/200. Solutions for kL/n in equation (fT5|) occur at intersections of 
the curves with a straight line through the origin with slope An/L. Thus the figure shows 
that in the limit |A|/L — )■ the degeneracy between the ordinary state and the lowest nodal 
energy at kL/ir = 5 is lifted by moving the potential. It also shows that, with i/L = 2/5 for 
example, the double degeneracy at kL/n = 5 (that occurs in the limit \A\/L — )■ 0) is split 
when A/L is nonzero and of either sign. The zeroes of the functions shown in Fig. [T] can 
easily be computed from Eq. (^^. 

To explore the rational case i/L = 2/5 a bit further we have calculated the lowest energy 
eigenfunction that is doubly degenerate in the limit \A\/L — )■ 0. The lowest energy nodal 
state is, from ( IT3l) . given by ^5" (x) = a/2/L sin(57rx/L). To get the ordinary solution at 
this energy we let kL = 5n + e in Eq. ([9]) and take the limit to find 

Dsm{5nx/L) for (2/5)L < x < L 
-(3/2)Dsin(57rx/L) for < x < (2/5)L, 



where D = 2/(v 3L) is a normalization constant. In the limit this ordinary solution vanishes 
at the delta function (x = 2L/5), suffers a discontinuity of slope there, and is orthogonal to 
the nodal solution ^5" . 

VI. THE SPECIAL CASE i = 1/2 

The simplest choice is to locate the delta potential precisely in the middle so that i/L is 
the rational number 1/2 and /i = 0. For this symmetric case it turns out that the ordinary 
and nodal states are precisely equal in number in the sense that their energies interleave. 
That is not the case for all other choices of i/L, where the symmetry is lost, but it is still 
true that both types are countably infinite in number and that their inclusion is necessary 
to make up a complete set. 

The nodal wave functions and their energies, '4'2ji^) ^^^ ^2,j, are given by ( IT3|) with 
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wave vectors kL = 2j vr, j = 1,2,3..., Eq. (fT^ . These energies i?2j are, of course, all 
positive. 

For the ordinary solutions any negative energy eigenvalue would have k = in and E = 
—fi^H? /{2m) where, putting yU = in fITTI) . k would satisfy the equation (kL/2)(2A/L) = 
tanh.{K,L/2). A sketch of this equation with respect to the variable x = nL/2 shows that there 
are no solutions for n when the interaction is a repulsion, none for relatively weak attraction 
(1 < 2A/L), and but a single solution for relatively strong attraction (0 < 2A/L < 1). In 
particular in the limit of very strong binding, for which 2A/L approaches zero from above, 
K —7- 1/A, and the energy approaches —Eb- 

For the ordinary solutions at positive energies, acceptable values for kL must, from (TTB]1 . 
obey the equation {kL/2){2A/L) = ta.n{kL/2). This will include all the ordinary solutions 
except when the potential is relatively strongly attractive, for which < 2A/L < 1 and the 
lowest energy moves to negative values. 

When the interaction is repulsive or relatively weakly attractive all ordinary states have 
positive energies, and a sketch reveals that, for this symmetric case, these solutions for kL 
interleave between those for the nodal solutions. In this sense, for the special case i = 1/2 
the nodal and ordinary states are equal in number. 

For a strong attraction (0 < 2A/L < 1) there is but one negative energy state (TT8l) . If 
the potential is repulsive, whatever its strength, then all states have positive energies. In 
addition to the nodal states ^2ji^) ^^ have the ordinary states. Of these, should there be 
attraction, there will be, as A/L -^ 0+, the single negative energy state with energy —Eb 
and wave function approaching \I/l/2, where (see Eq. (jH])) 

1 / |x-L/2r 
^^""l A 

For strong repulsion (and for strong attraction apart from the ground state "^1/2) all 
states have positive energies E = h k"^ /{2m). For these states \A\/ L is small and solutions 
for kL are in the neighborhood of the zeroes of tan {kL/2), namely kL = 2jtt + e, where e 
approaches zero as \A\/L — )■ and j = 1, 2, 3 . . .. Thus, as \A\/L — )■ 0, the energies of these 
ordinary states converge towards a double degeneracy with each and every one of the the 
nodal state values E2J, Eq. ( 1T3|) . except, for attraction, the lowest state. 

It is easy to find the positive energy wave functions for strong interaction. From flT3|) the 
nodal states are ^"2" (^) = \/2/L sin {2jttx/L). Using kL = 2jn + £ in Eq. ([9]), taking the 
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^L/2{x) = —j=. exp ( ^ ) for < X < L. 



limit £ — )■ and normalizing, gives the following expression for the ordinary states in the 
limit of strong interaction: 

^.(x) = e{x - (L/2))4"](x) (j = 1,2,3. . .), 

where 9 is the step function 

{ + 1 M >0 
(25) 
-1 M<0. 

These ordinary states have a discontinuous slope at the delta function, are even valued with 
respect to reflection about x = L/2 and are manifestly orthogonal to the nodal states which 
are odd valued about the middle. In this limit the energy of each ordinary state ipj{x) 
approaches that of the corresponding nodal state ipl^Ax): in this case of high symmetry all 
states approach doubly degeneracy in the limit of strong interaction. However we note that 
for all finite strengths of the contact potential all energies of this one- dimensional system 
are, of course, non-degenerate. 

VII. GROUND STATE 

We considered in subsection |Vl] solutions to Eqs. (ITB]) and ( 1T7|) with /x = 0, i.e. the 
case i = L/2. In that instance the ground state energy is positive for all strengths of 
repulsion, A < 0, and also for relatively weak attraction, 2A/L > 1. But for relatively 
strong binding, for which the inequality 1 > 2A/L > holds, the ground state energy is 
negative. And, in particular, in the limit of strong binding, when A/L — ?■ 0"*", it approaches 
— ^^/(2mA^) = —Eb, the bound state energy of the free system, Eq. (jS]). 

When the force center is not in the middle, then /i does not vanish and a re-analysis is 
required. For that case, whatever the details may be, we can say from the start that for 
either sign and any magnitude of A, because Eqs. (fT6!) and (TT7|) are symmetric with respect 
to fi about /i = 0, as a function oi i, any ordinary state energy eigenvalue must be symmetric 
about i = L/2. It is easy to see from f ll3p that this symmetry also applies to the nodal 
states. 

Consider, in particular, configurations with the delta function located near to the left 
wall, say, with i = e, where e is small and positive. As e decreases, the nodal solution 
energy eigenvalues — see Eq. (IT3l) — will rapidly increase to large positive values. To see 
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this, note that with i/L = p/n, n must increase without bound as e — )■ O"*", taking with it 
the accompanying energies Enj ~ (j nY, j = 1, 2, 3, ■ ■ ■ . 

As for the ordinary solutions, with i = e and e — )■ 0~^, we can say that as the delta 
function approaches a wall there are no negative energies. To see this, expand Eq. f lT7|) for 
small e to get tanh.{K,L) / (nL) ^ — 2e^/(AL). In the limit e — )■ 0^ this has no solutions for 
any finite value of k. So all energy eigenvalues are non-negative when the delta function 
potential is near to either wall. 

For arbitrary yU G (—1, 1), a sketch allows us to predict the existence of ordinary solutions 
for small positive and negative energies, the transition between which occurs at A; = 0. The 
second of Eqs. (12 2 p (which applies when kL is small) and the graphical picture suggests the 
following: that for small kL and an attractive potential (A > 0) there is a solution for small 
real valued kL provided the slope A/L is greater than (1 — yU^)/2; that k vanishes when 
A/L = (1 — /U^)/2; that there is a negative energy solution with k = in when A/L is less 
than (1 — /U^)/2. From Eq. ( fT6l) fi = 2i/L — 1, so we can say that as a function of i/L, the 
ground state energy passes through zero when A/L = [1 — /i^]/2 = 2(1 — i/L)i/L, or 

e/L = (1/2) (l ± y/1 - 2A/l) . (26) 

Note that, because i/L is restricted to the interval (0, 1), there can only be zeroes of energy 
provided < 2 A/L < 1, i.e. for relatively strong attraction. 

Finally, consider the ordinary solutions for values of i near, say, the left wall so that 
i = e, where e is small. Expanding Eq. ( !T5l) to lowest order and collecting terms gives 
tan{kL)/{kL) ~ -2t^/AL, so that in the limit t ^ Q, kL ^ Nix {N = 1,2,---), or 
E ^ h N'^Tx'^ /{2mL'^). Thus, in this limit the energies approach those of a free particle in a 
well, Eq. ([2]). To see how this limit is approached as £ — )■ 0^ we can set kL = A^7r + A, where 
A is small. Then tan{kL)/kL = tan A/(iV7r + A) ~ A/(iV7r), so that A/(iV7r) ~ -2eV(AL) 
or A ^ — 2A^7re^/(AL). The corresponding energy is 

This result holds for either sign of A = h'^ /(mX). The conclusion, then, is that when the 
delta function is near either wall, as a function of distance away, the energies are parabolic in 
form with zero slope at the wall, and they are a local maximum (minimum) for an attractive 
(repulsive) interaction. 
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The lowest, ground state, energy corresponds to the choice A^ = 1 in Eq. (127|) . Numerical 
calculation gives the dependence of energy eigenvalues on I as it is varied from to L. All 
energies above the ground state are positive so for them one must solve either of the pair f llSp 
or f ITB]) . For the ground state, when 2A/L > 1, these equations pertain, but when 2A/L < 1 
Eq. (I17p holds. Their solutions join smoothly at zero energy, for which k vanishes. This 
happens at values of C given by Eq. ( 126|1 . Results for the ground state are shown in Fig. [2] 
where the ratio E/Eb versus R = i/L is plotted for the three values (0.1,0.4,0.5) for the 
parameter / = A/L. We have done this for attractive and repulsive delta function potentials 
having the same strength |A|. For the former (latter) case A > (A < 0) and the the energy 
eigenvalues have their minimum (maximum) value at / = 1/2. 

Generally, for the case of attraction, should the wall separation L increase to values large 
with respect to the characteristic length A one expects that the ground state energy would 
approach —Eb, Eq. (|8]). Indeed in the limit kL -^ oo Eq. ( [T7|) becomes k -^ 1, provided 
< i < L so that K — )■ 1/A and E —^ —h /(2mA'^) = —Eb- Note, however, that however 
large L may be, according to Eq. (1271) with ra = 1, in the limits £ — t- 0^ and i — ?► L^ the wall 
dominates so that the energy approaches h 7r^/(2mL^), the lowest energy of a free particle in 
the well. This is consistent with Fig. [2l In the other extreme of small L, the walls dominate, 
whatever the value of i. In this limit |A|/L is large, which is equivalent to weak coupling. 
For this limit the analysis around the first of Eqs. f l22|) shows that kL ^ it as |A|/L — )■ cxd, 
so that in the limit L — t- we find that E ^ h n^ / (2mL^) = ir^f^EB, the lowest energy of 
a free particle in a well. 

VIII. CONCLUDING REMARKS 

The model analyzed here is simple and exactly solvable, but we are unaware of any com- 
plete published analysis. If we think of an attractive delta function as a very simple model 
of the nucleus of a relatively massive atom, then along the lines of the Born-Oppenheimer 
approximation we might suppose that its position i G (0, L) can be varied as a parameter. 
Thus the model suggests, reasonably, that for an attractive interaction the configuration 
of lowest energy is for the 'nucleus' to avoid the walls. This agrees with the perturbation 
calculations of Changa et al^ for hydrogen in a hard spherical cavity. Note that as the 
well width is decreased the uncertainty principle demands that eventually the kinetic energy 
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term of the Schrodinger equation is the dominant contribution to the energy. As we pointed 
out in Sec. I VIII this also applies for a well of any width should the fixed potential be located 
sufficiently close to either wall. 

We have shown that the energy levels and wave functions for choice i. = L/2 follow as the 
limit of a 'top hat' model potential having height H and width W < L where the 'area' HW 
is held constant and VT — )■ 0"*", and it is reasonable to conjecture that the same result applies 
for off-center top hat potentials and, more generally, to any sequence of model potentials 
which approaches the delta function as a limit. 

As a limit of sequences of non-singular potentials, the choice v{x) = —X6{x — i) does 
require care in its handling. In particular the behavior of the nodal energies, Eqs. (1131) and 
flT^ . is somewhat chaotic as i is varied. For instance, for any given 'rational' position i/L = 
p/n, the lowest nodal energy is proportional to n"^, so that the lowest such energies occur, 
in increasing order, at relative position i/L = 1/2 and then at i/L = (1/3, 1/4, 1/5, ■ ■ ■) 
to the left of the midpoint or, symmetrically, at points i/L = (2/3, 3/4,4/5, ■■ ■) to its 
right. But should i/L be moved through nodal points arbitrarily close to any one of these 
points then the lowest associated energy sweeps through chaotically large, even infinite, 
values. Nevertheless, we emphasize that, for rational values of i/L a complete set of solutions 
requires both the nodal and ordinary ones. By contrast, the ordinary solutions have energies 
which are solutions to Eqs. flTSl) or flT6|) . They are smooth continuous functions of position 
i G (0, L). For instance this is true for the ground state discussed in Sect. IVIII 

To conclude the paper we should like to make a few remarks of a purely mathematical 

nature concerning the model, most of which remains to be done. First we prove that the 

Hamiltonian, properly defined mathematically, is self adjoint. As the reader will know, self 

adjointness is that subtle property of a symmetric (hermitian) operator that guarantees that 

the operator can be exponentiated. For the Hamiltonian, this condition is necessary and 

sufficient to guarantee a unique time translation group. For our formal Hamiltonian, which 

we here denote 

S = -^^,-X6{x-i), 
2m dx'^ 

the process begins by defining a domain of definition for its interpretation as a sesquilinear 
form. We choose as domain T> the set of functions / G -^^(0, L) in the Hilbert space which 
belong to the domain of the kinetic operator and, in addition, have a first derivative /' which 
satisfies the discontinuity boundary condition ai x = i. This implies that /' G lv^(0, L). We 
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can then consider the form 



'9,f]--= 9it)[Sf]it)dt 

'0 



for all / in the domain V. By a simple calculation, involving no more than integration by 
parts, it follows that the associated quadratic form is lower bounded, 

I/,/l = ^ll/'ir-A|/Wr>-A sup |/(x)r>-oo. (28) 

'i'rn o<x<L 

It is now standard that the form can be extended to a closed form which defines a self 
adjoint operator, the Friedrichs extension of S. It is this extension that we take as the 
mathematically proper Hamiltonian H. For details of the extension we recommend Davies^^ 
and §124 of Riesz and Sz.-Nagy2^. 

This leaves the following to prove: the spectrum of H consists only of eigenvalues, indeed 
of those eigenvalues we have obtained earlier and only those. We firmly believe these to be 
true, and that, in consequence, for a given i and A the eigenf unctions we obtained constitute 
an orthonormal basis for L^(0, L). We remark that were we to independently prove that the 
eigenf unctions were complete then the eigenvalues we found would comprise the spectrum 
oiH. 
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Figures 
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FIG. 1: This shows the right-hand side (RHS) of Eq. (J15p as a function of kL in units of vr for 
i/L = 2/5 (dashed) and (./L = 0.415 (sohd). Solutions for the ordinary states are intersections of 
the curves with a straight hue from the origin having slope vrA/L. When (./L = 2/5 the lowest 
nodal state has a wave number given by kLj-K = 5, which is shared with an ordinary state (dashed 
line) in the limit |A| — t- 0. That the dashed curve moves from the value 5 shows that this limiting 
degeneracy is lifted as |A| moves from zero. On the other hand, note that even when A = 
this double degeneracy is lifted by moving the potential away from the value 2/5. The fact that 
0.415 = 83/200 is also rational does not effect this basic splitting: Its lowest lying nodal state 
occurs well to the right, at kL/n = 200. 
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FIG. 2: The result of numerical calculations (for values (0.1,0.4,0.5) of the parameter / = A/L) 
for the dependence of the 'Relative Energy' E/Eb of the lowest energy eigenvalue versus p = i/L 
locating the delta function within a box of side L. The curves with minima are for an attractive 
potential (A > 0). Their partner curves, with maxima, correspond to a repulsive interaction of the 
same strength. The energy of the single bound state of an attractive potential of strength |A| is 
E = —Ebi Eq. dH]), and A is the spatial decay length for that state, Eq. ([7]). 
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